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Abstract A classical and a rclativistic law of motion 
for a supernova remnant (SNR) arc deduced assuming 
an inverse power law behavior for the density of the in- 
terstellar medium and applying the thin layer approx- 
imation. A third equation of motion is found in the 
framework of rclativistic hydrodynamics with pressure, 
applying momentum conservation. These new formu- 
las are calibrated against a decade of observations of 
SN 1993J . The existing knowledge of the diffusive pro- 
cesses of ultrarelativistic electrons is reviewed in order 
to explain the behavior of the 'U' shaped profile of in- 
tensity versus distance from the center of SN 1993J. 

Keywords supernovae: general supernovae: individ- 
ual (SN 1993J ) ISM : supernova remnants 



1 Introduction 

The study of th e supernova remnant (SNR) started 
with lOorti (119461 ) where an on ongoing collisional ex- 
citation as a result of a post-explosion expansion of the 
SNR against the ambient medium was suggested. The 
next six decades where dedicated to the deduction of an 
analytical or numerical law of expansion. The target is 
a relationship for the instantaneous radius of expansion, 
R, of the type cx where t is time and m is a parame- 
ter that depends on the chosen model. On adopting this 
poi nt of v ie w, th e Sedov expansion predicts R cx t^'^, 
see Sedov (Il959l) . and the thin layer approximation 
in the presen ce of a constant density medium predicts 



-P- ^ 

R cx see lDvson. J. E. and Wilhams. D. A.I (|l997t) . 
A simple approach to the SNR evolution in the first 10^ 

L. Zaninetti 

Dipartimento di Fisica Generale, 
Universita degli Studi di Torino 
Via Pietro Giuria 1, 
1-10125 Torino, Italy 



yr assumes an initial free expansion in which R cc t un- 
til the surrounding mass is of the order of 1 Mq and 
a second phase characterized by the energy conserva- 
tio n in which according t o the Sedov solution R cx t^^^, 



see 



McCrav and Lavzerl ([1987|). A third phase charac- 



terized by an adi abatic expansion wit h R cx ^^Z*" starts 



after 10 yr, see iMcCrav and Layzeii (119871). A more 



sophisticated approach given bv lChevaliej ()l982lf 



alyzes self-similar solutions with varying inverse power 
law exponents for the density profile of the advanc- 
ing matter, R~", and ambient medium, R^^. The 
previous assumptions give a law of motion R cx t"-= 
when n > 5. An other example is an a nalyt ical so- 
lution suggested by Truelove and McKee ( 19991) where 
the radius-time relationship is regulated by the de- 
crease in density: as an example, a density proportional 
to R^^ gives R cx t^^^. With regard to observations, the 
radius-time relationship was clarified when a decade of 
very-long-baseline interferometry (VLBI) observations 
of SN 1993 J at wavelengths of 3.6, 6, and 18 cm became 



available, see iMarcaide et al\ ( 20091 ). As a first exam 



pie, these observations collected over a 10 year period 
can be approximated by a power law dependence of the 
type R cx This observational fact rules out the 

Sedov model and the momentum conservation model. 
The observed radius-time relationship leaves a series of 
questions unanswered or merely partially answered: 

• Is it possible to deduce a classical equation of motion 
for the SNR with an adjustable parameter that can 
be found from a numerical analysis of the radius-time 
relationship? 

• Is it possible to deduce a relativistic-mechanics equa- 
tion of motion for the SNR, since the initial velocity 
of the SNR can be on the order of 1/3 of the velocity 
of light? 

• Is it possible to deduce a relativistic-hydrodynamics 
equation of motion for the SNR applying momentum 
conservation? 
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• Can we build a diffusive model which explains the 
behavior of the emission intensity of the SNR? 

• Can a simple model for the time evolution of the total 
mapped flux densities of SN 1993J be built? 

In order to answer these questions, Section [2] reports 
the data on SN 1993J . Section [331 reports a new classi- 
cal law of motion assuming an inverse power law depen- 
dence for the density of the medium. Section 2] reports 
the evolution of SN 1993J on the basis of four models. 
Section [5] contains two new relativistic laws of motion 
as well as a fit to the data. Section [6] reviews the exist- 
ing situation with the radiative transport equation as 
well as three different processes of diffusion for ultrarel- 
ativistic electrons. Section |6] also contains a new scaling 
law for the temporal evolution of the flux of SN 1993J at 
6 cm. 



2 A spherical SNR 

The supernov a SN 1993J st a rted t o be visible in M81 

(jl993l ). a nd presented a 



in 1993, see iRipero et al 



circula r symmetry for 4000 days, see iMarcaide et al. 



( 20091). Its distance is 3.63 Mpc (the same as M81), see 



Freedman et al\ (|l99-f ). The expansion of SN 1993 J has 



been monitored in various bands over a decade and Fig. 
[1] reports its temporal evolution. 

The instantaneous velocity of expansion can be de- 
duced from the following formula 



tj+i - tj 



(1) 



where rj and tj denote the radius and the time at 
the position j. The uncertainty in the instantaneous 
velocity is found by implementing the error prop- 
agation equation (often called the law of errors of 
Gauss) when the covariant terms are negle cted, see 
Bevington. P. R. and Robinson. D. kI |2003f )). Fig. [2] 
reports the instantaneous velocity aw well the relative 
uncertainty. In particular, the observed instantaneous 
velocity decreases from v = 15437 ^ at i = 0.052 yr 



3.1 The constant expansion velocity 

The SNR expands at a constant velocity until the sur- 
rounding mass is of the order of the solar mass. This 
time , tj\/ , is 



tM = 186.45- 



'^0 'VlOOOO 



yr 



(2) 



where Mq is the number of solar masses in the volume 
occupied by the SNR, no, the number density expressed 
in particles cm""^, and Wioooo th e initial velocity ex- 
presse d in units of 10000 km/s , see lMcCrav and Lavzer 
(Il987t) . 



3.2 Two solutions 

A first law of motion for the SNR is the Sedov solution 
'25 



i?(t) 



TT p 



(3) 



where E is the en e rgy in j ected into the process and t 
is time, see Sedov ( 19591) : McCrav and Lavzer (1987). 
Our astrophysical units are: time, (ti), which is ex- 
pressed in years; -E51, the energy in 10^^ erg; uq, 
the number density expressed in particles cm~'^ (den- 
sity p = nom, where m = 1.4toh)- In these units, equa- 
tion ([3]) becomes 



R{t) « 0.313 



ti' 



no 



pc 



(4) 



The Sedov solution scales as t" '*. We are now ready to 
couple the Sedov phase with the free expansion phase 



0.0157 tpc 



0.0273 V^Vc 



if t < 2.5yr 
if t > 2.5yr 



(5) 



This two phases solution is obtained with the following 
parameters A/0 =1 , no = 1.12710^ , -E51 = 0.567 and 
Fig. |3] reports it's temporal behavior as well the data. 
The quality of the simulation , e„, can be obtained by 



to V = 8474 ^ at t = 10.53 yr. We briefly recall that a comparison of observed and simulated quantities: 



Fransson et al\ (j2005l ) quote an inner velocity from the 



shapes of the lines of 
of « 10000^. 



3 Classical case 



7000^ 



and an outer velocity = (1 — 



|(-Rpc,obs Rpc.num)\ ' 



R 



100. 



pc.obs 



(6) 



where Rpc,oha is the observed radius, in parsec and 



Rr, 



is the radius from our simulation in parsec. In 



This section reviews the free expansion , two simple 
laws of motion for the SNR and a new law of motion in 
the light of the classical physics. 



the case of the two-phases simulation we have tq = 65% 
which is a low value in comparison with the other mod- 
els here considered see Table [TJ 

A second solution is connected with momentum 
conservation in the presence of a constant density 
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Time (yr) 



Fig. 1 Radius in pc versus year of SN 1993J witli vertical error bars. Tiie em bands are A = 3.6 cm and A — 6 cm. The 
data are extracted from Table 1 in Marcaide et al. 2009. The dotted line represents an expansion at a constant velocity. 




Time 

Fig. 2 Instantaneous velocity of SN 1993J with uncertainty. 
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Time (yr) 

Fig. 3 Theoretical radius as given by the two-phases solution with data as in Table [T] (full line) , and astronomical data 
of SN 1993J with vertical error bars. 
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medium, see Dyson. J. E. and W il liams. D. A.I (|l997l ): 
PadmanabhanI |200l[) : IZaninettil (|2009l) . The astro- 
physical radius in pc as a function of time is 



R{t)= {jRo''{AmiQ~^vi {h-to) + Ro)pc , (7) 
where ti and <o arc times in years, i?o is the radius in 



1 fer; 



units 



pc when ti = to and vi is the velocity in 
when ti = to. The momentum solution in the presence 
of a constant density medium scales as t° '^^. 

3.3 Momentum conservation with variable density 

We assume that around the SNR the density of the in- 
terstellar medium (ISM) has the following two piecewise 
dependencies 



po if i? < Ro 

poij^)'' iiR>Ro 



(8) 



In this framework, the density decreases as an inverse 
power law with an exponent d that can be fixed from the 
observed temporal evolution of the radius, with d = 
meaning constant radius. The mass swept, Mq, in the 
interval < r < i?o is 



Mo = - PoTT ^0 



(9) 



The mass swept, M, in the interval < r < R with 

> -Ro is 



M 



PottRo^ 4 3 



Momentum conservation requires that 
Mv = MoVo , 



(10) 



(11) 



where v is the velocity at t and vq is the velocity at 
t = to- This formula is not invariant under Lorentz 
transformations and the initial velocity, vo, can be 
greater than the velocity of light. The previous ex- 
pression as a function of the radius is 
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rp^Po (3-d) 



(12) 



where l3o=vo/ci, P=v/ci and q is the velocity of light. 
Here, we have introduced a relativistic notation for later 
use. In this differential equation of first order in R, the 



variables can be separated and an integration tcrm-by- 
term gives the following nonlinear equation J-nl 

Tnl = (4i?o^rf - Ro^d^) R - iRo'^R^^'^ + Ro^d^ 
+12Ro^vot + 3i?o* - 4i?o*d 
+7Ro^vodto + Ro^vod^t 
-IRo^vodt - URo^votg - Ro^vod^to = . (13) 

An approximate solution of J-nl{t) can be obtained 
assuming that ZR'^R'^-'^ > -(4i?i^rf - Rf^d'^)R 



(i?o 



4-(i 



dRo^-'^ii - d) 



+ -(4 - d)i>oi?o'-''(3 - d)it - to)) — 



(14) 



Up to now, the physical units have not been specified, 
pc for length and yr for time are perhaps an accept- 
able choice. With these units, the initial velocity vq 
is expressed in — and should be converted into — ; 
this means that vo = 1.02 10~^wi where vi is the initial 
velocity expressed in 

The astrophysical version of the above equation in 
pc is 



R{t) = 

{Ro^-'^ - \dRo^-'^{A -d) + 3.402 10~^ x 
x(4-d)wii?o'"''(3-d)(ti-io))^pc , 



(15) 



where ti and to are times in years, Rq is the radius in 
pc at ti = to and vi is the velocity at ti = to in — . 



4 Classical fits 



The quality of the fits is measured by the merit function 



(Rth — RobsY 



(16) 



where Rth, Robs and aobs are the theoretical radius, 
the observed radius and the observed uncertainty re- 
spectively. 

A first numerical analysis of the observed radius- 
time relationship of SN 1993J can be done by assuming 
a power law dependence of the type 



R{t) = rpt"" 



(17) 



The two parameters and ap can be found from the 
following logarithmic transformation 



MR{t))=\n{rp)+a.p\n{t) , 



(18) 
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which can be written as 



y = a-LS + hhsx 



(19) 



The apphcation of the least square niethod through the 
FORTRAN subroutine LFIT from IPress et al\ (l992l) 
allows to find a^g Jd^s a-nd the errors a a and ab, 
sec numerical values in Table [1] The error on Tp is 
found by implementing the error propagation equa- 
tion (often called law of errors of Gauss) when the 
co yariant terms are neglected (see eq u atioii (3.14) 



Bevington. P. R. and Robinson. D. k1 |2003l) ). 



(Tr-p = exp(a)cro 



(20) 



and (7ctp = cTfe. In this case, the velocity is 
V^(t)=rpapi("--i) . (21) 
A second numerical analysis can be done by assum- 



ing a piecewise function as in Fig. 4 of iMarcaide et 
(I2OO9I) 



R{t) = 



(22) 



This type of fit requires the determination of four pa- 
rameters, i.e., tbr, the break time, r^r the radius of ex- 
pansion at t = thr and the exponents of the two phases 
are alphai and alpha2. The two- regime fit can be visu- 
alized in Fig. Uland Table[T]reports the four parameters. 

A third type of fit can be done by adopting the ap- 
proximate radius as given by equation (jisp ; Fig. [5] re- 
ports a fit of this type. 

A fourth type of fit implements the nonlinear equa- 
tion ([13]); the four roots can be fo und with th e FOR- 
TRAN subroutine ZRIDDR from IPress et al\ (|l992[) : 
Fig. ini reports a fit of this type. 



5 Relativistic case 

The relativistic analysis is split in two. 

1. The relativistic mechanics case in which the temper- 
ature effects are ignored and the thin layer approxi- 
mation is used. 

2. The case of relativistic hydrodynamics in which the 
pressure is considered and the velocity is found from 
momentum conservation. 



Table 1 Numerical values of the parameters of the fits 
and x^- ^ represents the number of free parameters and €q 
the quality of tlie simulation. 



N 



vah 



power law 



ap ^ 0.82 ± 0.0048 



(0.015 ±0.00011) 



6364 98.54% 



piecewise 

ai = 0.83 ±0.01 
^2 = 0.78 ±0.0077; 

rbr = 0.05 pc; tbr — 4.10 yr 



32 97.76% 



approximate radius 
rf = 2.54;ro = 0.019 pc; 
to = 0.249 yr;vo = 100000-^ 



7186 96.95% 



nonlinear radius 
4 rf = 2.93; ro = 0.019 pc; 276 93.3% 
to = 0.249 yr;vo = 100000-^fi 
relativistic radius 
4 d = 2.54; ro = 0.0045 pc; 
to = 0.249 yr;l3o = 0.333 



5557 93.05% 



relativistic pressure radius 
5 d = 0.89; ro = 0.00072 pc 

/ = 10-5 

to = 0.249 yr-l3o = 0.333 



1844 99.93% 



5.1 Relativistic mechanics 



Newton's law in special relativity, after lEinsteinl (|1905[ ). 
is: 



(23) 
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2 4 6 a 10 



Time (yr) 

Fig. 4 Theoretical radius as given by the two- regime fit represented by equation (|22p with data as in Table [T] (full line), 
and astronomical data of SN 1993J with vertical error bars. 
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Fig. 5 Theoretical radius as given by the approximate radius, equation l|15p. with data as in Table [T] (full line). The 
astronomical data of SN 1993J are represented with vertical error bars. 




Fig. 6 Theoretical radius as obtained by the solution of the nonlinear equation l|13p (full line), data as in Table [T] The 
astronomical data of SN 1993J are represented with vertical error bars. 
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with 



5.2 Relativistic hydrodynamics 



1 



(24) 



where F is the force, p is the relativistic momentum, m 
is the relativistic mass, is the rest-mass, v is the ve- 
lo city and c; is the vel ocity of light, see equation (7.16) 



French. A.P.I (Il968l) . In the case of the relativistic 
expansion of a shell in which all the swept material re- 
sides at two different points, denoted by radius R and 
radius i?o, the previous equation gives: 



M- 



Mn 



(25) 



where l3o=vo/ci, f3=v/ci, M is the rest mass swept be- 
tween and R and Mq is the rest mass swept between 
and i?o- This formula is invariant under Lorentz trans- 
formations and the initial velocity, vq , cannot be greater 
than the velocity of light. Assuming a spatial depen- 
dence of the ISM as given by formula ([U, relativistic 
conservation of momentum gives 



3 (-3 + d) y/T^ 3 Vl-l3o'' 

According to the previous equation, /3 is 

-i-3 + d) f3oro^ 



with D 



/3 



13 



(26) 



(27) 



-Qr^-'^Ra 



6-2d n 2d 

'^+^d + Ro^d^ 



-9/3o 



Ri 



2d 



6 /3o r^^^i^o d 



In this differential equation of first order in r, the vari- 
ables can be separated and the integration can be ex- 
pressed as 



-Ro 



^dR = c (3 - rf) (3oRo^ [t - to) . 



(28) 



The integral of the previous equation can be performed 
analytically only in the cases d = 0, d = 1 and d = 2, 
but we do not report the result because we are inter- 
ested in a variable value of d. The integral can be eas- 
ily evaluated from a theor etical point o f view using the 
subroutine QROMB from lPress all (|l992f) . The nu- 
merical result is reported in Fig. [7] and the input data 
in Table [TJ The behavior of relativistic and classical 
velocities are reported in Fig. HI 



A relativistic flow on flat space time is described by the 
energy-momentum tensor, T^'^ , 



pg 



(29) 



where is the 4-velocity, and the Greek index 
varies from to 3, w is the enthalpy for unit vol- 
ume, p is the pre ssure and g^'^ t he inverse met - 
ric of the manifold \yeinberg ( 19721): Landau ( 1987 ) 



( 20051) : 



Gourgoulhon 



2006 ) 



Hidalgo and Mendoza 
Momentum conservation in the presence of velocity, 
V, along the radial direction states that 



(z^(-) 



f s2 1 



Q 1 



p)A = cost 



(30) 



where A{r) is the considered surface area, which is per- 
pendicular to the motion. The enthalpy per unit vol- 
ume is 



ClP + P 



(31) 



where p is the density, and ci is the velocity of light. The 
reader may be puzzled by the 7^ factor in equation pop , 
where 7^ = "'"..a . However it should be remembered 



that w is not an enthalpy, but an enthalpy per unit 
volume: the extra 7 factor arises from 'length contrac - 
tion' in the direction of motion Gourgoulhon ( 20061 ). 
We continue assuming 



p : 



(32) 



where / is a parameter which has the range < / < 1 
and is supposed to be c onstant during t he expansion, 
see formula (2.10.26) in IWeinberd (|l972l ) for a hot ex- 
tremely relativistic gas. The previous equation (|30l) 
becomes 



A = constant . (33) 



The density is supposed to vary during the expansion 
as 




f^Osd 

P = Poi^) 



(34) 
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Fig. 7 Theoretical radius as obtained by the solution of the relativistic-mechanics equation (|28|) (full line), data as in 
Table [T] The astronomical data of SN 1993J are represented with vertical error bars. 




Fig. 8 Theoretical relativistic-mechanics velocity as given by equation (|27|) (dashed line), theoretical classical velocity as 
given by equation (|12|l (full line) and instantaneous velocity of SN 1993J with uncertainty. 
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where po is the density ai R = Rq. In two surfaces of 
the expansion we have: 



1/3 R'f{^l + I3T' 

1/3 



(35) 



The positive solution of the second degree equation is: 

N 



N 



D 

+ 6 

-3i?o"'i?"''+'/3oV 



(36) 



-6 i?o'+'i?'''+'/3o'/ + 9 Ro^+^R-^+^po^ 
-3 i?o' ""R-^ V + 3 Ro'^+^R-^+^f 
+9 i?o'/3o'' + 3 i?o''+'/3o''i?"'+'/ + 6 i?o^/3o'/ 
-i?o''+V'^"'+' + Ro'^+^fR-'^+^Po^ + Ro\f)-- 

D = 3Ro'^R^'^+^Po^ 
3i?o^/3o^ - Ro^f ■ 



-3 Rq^R 



d+2 



The equation of motion is 
/ -di? = c(t-io) • 



(37) 



The i ntegral is evaluate d using the subroutine QROMB 
from lPress et ali (|l992l ). The numerical results are re- 
ported in Figs 1^ and [TUl Fig. [TU] reports the decrease 
of the relativistic velocity. 



6 How the image is formed 

In this section, the existing knowledge about adiabatic 
and synchrotron losses, the acceleration of particles 
by the Fermi II mechanism, , 3D mathematical diffu- 
sions with constant diffusion coefficients and the rim 
model with constant density are reviewed and applied 
to SN 1993J . A new example of a ID random walk with 
a step length equal to the relativistic electron gyro- 
radius is also reported. A new simple model for the 
temporal evolution of the flux densities is introduced. 

6.1 Acceleration and losses 

6.1.1 Adiabatic losses 

An ultrarelativistic gas which experiences an expansion 
loses energy at the rate 



,dE, 1,„ 



(38) 



where E is the energy and V • v is the diver gence of 
the ex pansion velocity, see formula (11.27) in iLongaiii 
( I994I ) . A simple expression for V • v can be found from 
the power law model, sec equations ([T71) and (pi]) 



V- V 



R- 



(3a- 1) 



R 



(39) 



where R is the temporary radius of the expansion and 
Vp and a are reported in Table [TJ The lifetime. Tad, of 
an ultrarelativistic electron for adiabatic losses is 



Tad = ^ = 344.39 i?i-22yr 

dt 



(40) 



when the radius R is expressed in pc. During the ten 
years of observed expansion, the radius of SN 1993 J has 
grown from « 0.01 pc to « 0.1 pc and therefore the time 
scale of adiabatic losses has increased from « 1.25 yr 
to « 20.7 yr. 

6.1.2 Synchrotron losses 

An electron which loses its energy due to synchrotron 
radiation has a lifetime r^, where 



(41) 



E is the energy in ergs, H the magnetic field in Gauss, 
a nd Pr is the t otal radiated power, see formula (1.157) 
inlLangI |l999l ). 

The energy is c onnected to t he critical frequency, see 
formula (1.154) in lLangl |l999l ). as 



Vc = 6.266 X lO^^HE^ Hz . (42) 
The lifetime, Tgym for synchrotron losses is 
1 



Tsyn = 39660 



yr 



(43) 



The time-scale of synchrotron losses is shorter than that 
of the adiabatic losses if the following inequality is ver- 
ified 



1.572 10^ 



(44) 



where H is expressed in Gauss and Tad in years. The 
previous equation can also be expressed as 



H > 1162.98 



1 



^tad'/^ 



Gauss 



(45) 



that is the magnetic field at which the synchrotron 
losses prevails on the adiabatic losses and Fig. [TT] re- 
ports the numerical values of the transition. 
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QCI o k_i , I , I • I , I , |_ 

2 4 • a 10 



Time (yr) 

Fig. 9 Theoretical radius as obtained by tiie solution of the relativistic-hydrodynamics equation (|37|l (full line) and data 
as in Table [T] with uncertainty. 




R (pc) 



Fig. 10 Theoretical relativistic-hydrodynamics velocity with pressure as given by equation (|36p (full line), and instanta- 
neous velocity of SN 1993J with uncertainty. 



Magnetic field of transition (Gauss) 




20 1. y 10^" 



Fig. 11 Magnetic field in Gauss over which the synchrotron losses prevail on the adiabatic losses. 
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6.1.3 Particle 



deration 



Fermi 


( 


1949, 


1954) 



of a particle which spirals around a line of force is pro- 
portional to its energy, E, 



B^E 



(46) 



where B = u/ci , see formula (3) in iFermil (|1949I ). The 
continuous form is 



dE _ E 



(47) 



where r is the typical time-scale. The probability, P{t) 
, that the particle remains in the reservoir for a period 
greater than t is now introduced, 



Pit) 



e T 



(48) 



where T is the time of escape from the considered re- 
gion. The resulting probability density, N{E), is 



N{E)^^{^r-'f , 

where Eq is the initial energy and 

7/ = 1 

Equation (UHl) can be written as 

N{E) = K E-^f , 



r 
T 



(49) 



(50) 



(51) 



where K 



A power law spectrum in the 



particle energy has now been obtained. In Fermi II 
processes, the typical time-scale, r//, when the particle 
stays in the accelerating region a time greater than T 
is 



TJJ 6 cf L 



(52) 



where u is the velocity of the accelerating cloud and L 
is the mea n free path b etween clouds, see formula after 
(4.439) inlLangl (Il999l) . The mean free path between 
the accelerating clouds in the Fermi II mechanism can 
be found from the following inequality: 



Tll < Tsyr, 



or 



L < 1.72310^ 



pc 



(53) 



(54) 



As an example, inserting /? = 2.82 10 ^ (value of ve- 
locity at « 0.1 pc of SN1993J ), = 1.5GHz and 



H = 65.1 Gauss (Table 1 in iMarti-Vidal et all (|201C 
we obtain 



L < 6.7610" 



pc 



(55) 



This model of acceleration can work in the rim model 
with constant density of emitting particles, see Section 
16.41 In this case the thickness of the emitting region is 
b — a=0.035 pc which means that an high number of 
collisions can be done. In Fermi II process the energy 
increases exponentially with time 



E{t) = Eo cxp(i/T) , 



(56) 



see , equation (3) in iFermil (jl954l ). The effect of the 
synchrotron losses on the energy of the electron dur- 
ing the various collisions and the consequent equilib- 
rium energy has been analyzed in IZaninetti and Siah 
(jl988() . The strong shock acceleratin g me c hanism , 
named Fermi I, was later introduced by iBell (Il978l l?h 
and produces an increases in energy of the particle of 
the type 



AE = -E 
ci 



(57) 



w here u is the velo city of the shock, see formula (21.20) 



Longair ( 1994h . This process allows to predict a 



probability density in the energy of the accelerated par- 
ticles of the type N{E) oc E""^. In our case can be 
considered the process which accelerates the particles 
before they start to diffuse from the shock , see Section 
1 6.31 and 16.41 A modern review o f the two Ferm i mech - 
anisms can be found in jKulsrudI dioosh : ISomovl (l2006l) : 
Dermer and MenonI (|2009l ). 



6.2 The transfer equation 

The transfer equati on in the presence of emiss ion 
only, see for example Rvbicki and LightmanI (|l99l[ ) or 
Hiellming. R. M.I (|l988l) . is 



— — = -k^plu 
ds 



-JuP 



(58) 



where I^, is the specific intensity, s is the line of sight, 
ji, is the emission coefficient, is a mass absorption 
coefficient, p is the density of mass at position s and the 
index v denotes the frequency of emission of interest. 
The solution to equation ([55)) is 



/.(T.) = ^(l-e--''(^)) , 



(59) 



where is the optical depth at frequency v 

dr^, = kypds . (60) 
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The volume emissivity (power per unit frequency in- 
terval per unit volume per unit solid angle) of the ul- 
trarelativi stic radiat ion from a group of electrons, ac- 
cording to iLang IS 



e{iy) = / P{v)N{E)dE 



(61) 



where P{v) is the total power radiated per unit fre- 
quency interval by one electron and N{E)dE is the 
number of electrons per unit volume, per unit solid an- 
gle along the line of sight that are moving in the direc- 
tion of the observer and whose energies lie in the range 
E to E + dE. In the case of a power law spectrum. 



N{E)dE ^ KE' 



If 



(62) 



where i^T is a constant. The value of the constant K 
can be found by assuming that the probability density 
function f or the relativistic energy is of Pareto type as 
defined in Evans. Hastings, and PeacodJ ( 20001 ) 



(63) 



with Cp > 0. In our case, Cp = 7/ — 1 and a = Emin, 
where Emin is the minimum energy. We can now ex- 
tract 

K ^ N„i^f - 1)EZ-' , (64) 

where Nq is the total number of relativistic electrons per 
unit volume, here assumed to be approximately equal 
to the matter number density. The previous formula 
can also be expressed as 



K = 



P 



1.4 TTiH 



(65) 



where mu is the mass of hydrogen. The emissivity 
of the ultrarelativistic synchrotron radiation from a ho- 
mogeneous and isotropic distribution of electrons whose 
N(E) is given by equation ([62]) is, according to iLang 

mm, 



0.933 X lQ-^^a{-if)KH 

,6.26 X IQi^ 



r(7/-l'l)/2 



(66) 



(7/-l)/2 



erg sec ^cm ^Hz ^rad ^ , 

where v is the frequency and a (7/) is a slowly varying 
function of 7/ which is of the order of unity and is given 

by 



a (7/) 

of-yf-3)/2 7/ + 7/3 . 37/ ~ l Np/37/ + 7.^ 
7/ + 1 ^ 12 ' ^ 



12 



(67) 



for 



> 



7/ ^ 2- 

We now continue to analyze the case of an optically 
thin layer in which is very small (or k^, is very small) 
and where the the density p is substituted for the con- 
centration C(s) of relativistic electrons 



(68) 



where is a constant function of the energy power law 
index, magnetic field and frequency of e.m. emission. 
The intensity is now 

I^{s) = Kg f C{si)dsi Optically thin layer . (69) 

J So 

The increase in brightness is proportional to the con- 
centration integrated along the line of sight. In numeri- 
cal experiments, the concentration is memorized on the 
visitation-grid S and the intensity is 



k 

Optically thin layer , 



(70) 



where As is the spatial interval between the various 
values and the sum is performed over the interval of 
existence of index k. The theoretical flux density is 
then obtained by integrating the intensity at a given 
frequency over the solid angle of the source. In order 
to deal with the transition to the optically thick case, 
the intensity is given by 

/(*,j) - ^(l-exp(-/^,^ As X S[i,j,k))) (71) 



k 

Thin 



Thick 



where Ka is a constant that represents the absorption. 
Considering the Taylor expansion of the last formula 
((7T|) . equation ([70| is obtained. 

6.3 3D diffusion from a spherical source 

Once the concentration, C . and diffusion coefficient, D. 
are introduced, Fick' s law in three dimensions is 



dC_ 



Under steady-state conditions, 
DV^C = . 



(72) 



(73) 



The concentration rises from at r=a to a maximum 
value Cm at r=h and then falls again to at r=c. The 
solution of equation ([73]) is 



C(t)=A+- , 
r 



(74) 
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where A and B are deterniined by the boundary condi- 
tions, 



The third segment, I^^\y), is 



a,(r) = C™(l-^)(l-^)"' a<r<h , (75) 6C„cln(y2) 5a„ch, ( + c) 



and 



cur) - a„ - 1) - 1) ' ^ 



(76) 



These solutions can be found in lBerg|(ll993l) or in lCranlj 
(|l979l ). Fig. Oshows a spherical shell source of radius h 
between a spherical absorber of radius a and a spherical 
absorber of radius c. The concentration rises from at 
r=a to a maximum value C,„ at r=b and then falls 
again to at r=c. The concentrations to be used are 
formulas (|75p and (j76p once r = \J -\- is imposed; 
these two concentrations arc inserted in formula (|59p 
which represents the transfer equation. The geometry 
of the phenomenon fixes thre e different zones (0 — a, a — 
b,b — c) for the variable y, see Zaninetti ( 2007t ): the first 
segment, I^{y), is 

I'iy) = 

^bC^^a-^-y^ ^bCraa\n(^^a2 -y^ + 



-b + a 



-b + a 



_„&WEZZ ^ ^ ^^rna\n[^b^-y^ + b) 
—0 + a —0 + a 



-c + 6 -c + 6 



— 



< y < a 



(77) 



The second segment, I^^{y), is 



^^^^^^ ^ bCrr^aln {y^) _ ^bCn^b^-y 



—b + a —b + a 

&C™aln(y6^^ + 6) 



+2 



-6 + a 



-c + b 



2 5 CmC In 

- 2 



— y^ + c 



-c + b 

-c + b 
a < y < b . 



(78) 



-c + 6 



-c + b 



+2 



-C + b 
b < y < c 



(79) 



The profile of / made up of the three segments ([77)) . 
(|75|) and ([7^ . can be calibrated against the real data 
of SN 1993J and an acceptable match can be achieved 
by adopting the parameters reported in Table [H The 
theoretical intensity can therefore be plotted as a func- 
tion of the distance from the center, see Fig. [T31 or as 
a contour map, see Fig. [Ml 

The position of the minimum of / is at y = and 
the position of the maximum is situated in the region 
a < y < b, or more precisely at: 



y 



^ — {b — 2a + c) a [ab — 2bc + ac) 
b-2a+c 



(80) 



This means that the maximum emission is not at the 
position of the shock, identified here as 5, but shifted a 
little towards the center; see Fig. [T51 

The ratio between the theoretical maximum inten- 
sity, Imax, as given by formula (|80p . and minimum in- 
tensity {y = 0) is complex and is reported in formulas 
(74-76) in lZaninettl (|2007l ). The observed ratio as weh 
as the theoretical ratio are reported in Table [2] 

The effect of absorption is easily evaluated by apply- 
ing formula ([7T|) and fixing the value of Ka- The result 
is shown in Fig. 1161 



Table 2 Simulation of SN 1993J by 3D diffusion, optically 
thin case 



symbol meaning 



value 



a radius internal sphere 1.76(mas) 

b radius shock 2.2{mas) 

c radius external sphere b.Q{mas) 

/''"*'' ratio observed intensities 1.7926 

ratio theoretical intensities 1.7927 
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Fig. 12 The spherical source is represented by the dashed line and the two absorbing boundaries by full lines. The 
observer is situated along the x direction, and three lines of sight are indicated. Adapted from Fig. 3.1 in Berg (1993). 



o 
E 



Radius (mas) 

Fig. 13 Cross-section of the mathematical intensity / (formulas (|77p . (|78|l and (|79|l 'l. through the center (dotted line) of 
SN 1993J and real data (empty stars), = 100.49 The real data made on day 1889 after the explosion have been extracted 
by the author from Fig. 3 of Marcaide et al. (2009). Parameters as in Tabled 




I 



mas from SNR center 

Fig. 14 Contour map of / adjusted to simulate SN 1993 J . Parameters as in Table (2] 



Models for SN 1993J 



15 




1 2 3 4 B 

Radius (mas) 



Fig. 15 Cross-section tiirougii tlie matiiematical intensity / towards tlie edge of SN 1993J . Tlie tliree parameters wliicli 
ciiaracterize tlie expanding PN, a , b and c, are reported. Parameters as in Table [5] 




Fig. 16 Cross-section through the mathematical intensity / (formulas (|77p . (|78|1 and (|79p \ in the optically thin case 
(dashed line, = 237.3), and optically thick case (full line, ~ 84.7) and real data (empty stars). Parameters as in 
Tabled 
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6.4 The rim model with constant density 



and 



We assume that the number density of ultrarelativistic 
electrons C is constant and in particular rises from 
at r = a to a maximum value Cm , remains constant 
up to r = b and th en falls again to 0, see Section 5.2 
in IZaninettil (|2009l ). The length of sight , when the 
observer is situated at the infinity of the x-axis , is the 
locus parallel to the a;-axis which crosses the position y 
in a Cartesian x—y plane and terminates at the external 
circle of radius b. The locus length is 



U = 2 X ( - 2/2) ;a<y<b . (81) 

When the number density of ultrarelativistic electrons 
C'm is constant between two spheres of radius a and b 
the intensity of radiation is 

loa = 



C„,x2x{^U^-y^-^a^-y^) ;0 < y < a (82) 
Iab^CrnX2x {^b^ - y^) ; a < y < b . (83) 

The ratio between the theoretical intensity at the max- 
imum , (y = a) , and at the minimum , (y = 0) , is 
given by 



l{y ^ a) ^ VP 
I{y = 0) 



b — a 



(84) 



The parameter b is identified with the external radius 
of the SNR. The parameter a can be found from the 
following formula 



^i(y 


=0) 




O) X 


i{y= 


0) > 


i{y= 
rr,,- 


--a) \ 



obs 



1 



(85) 



is the observed ratio between maxi- 
mum intensity at the rim and intensity at the center. A 
cut in the theoretical intensity is reported in Figure [TTI 

6.5 ID diffusion 

The Fick eq uation i n ID with a constant diffusion co- 
efficient , see ICrankI (Il979l) . IS 







The general solution to equation ((86)) is 

C{r)^A + Br . 

The boundary conditions give 



C(r) = C,n^ a<r<b 



— a 



(86) 
(87) 
(88) 



C{r) = C„ 



b < r < c 



(89) 



The transport of relativistic electrons with a step length 
equal to the gyro-radius of the relativistic electrons 
is call ed Bohm diffusion ( Bohm. Burhop. and MassevI 
(Il949l) ) and the diffusion coefficient is energy-dependent.| 
The assumption of Bohm diffusion allows of setting a 
one-to-one correspondence between the energy and the 
step length in the random walk. The relativistic elec- 
tron gyro-radius rn is 



meClV±Je 

qB 



(90) 



where rrie is the electron mass, gammae is the Lorentz 
factor of the relativistic electron, q is the velocity of 
light, Vperp is the velocity perpendicular to the mag- 
netic field, B is the magnetic field in Gauss and q is 
the electron charge in statco ulomb s . The astrophysical 
version, see formula 1.153 in iLang ( 1999t ). is 



TH = 210^— cm 
B 



(91) 



where E is the electron energy in cgs. The typical fre- 
quency of synchrotron emission, j^c, is 



6.266 lO^^SS^ i?2 



(92) 



This formula allows us to express the relativistic elec- 
tron gyro-radius as a function of the wavelength of emis- 
sion expressed in centimeters, Ai, 



1.4210" 



-PC , (93) 
5, is expressed in units of 



where the magnetic field, iJ_ 
10-5 Gauss. 

The ID theoretical solution and a Monte Carlo simu- 
lation characterized by a given number of trials, NTRI- 
ALS, are reported in Fig. [181 The theoretical intensity 
can therefore be evaluated by a numerical integration 
of equations (|88|) and ([89]) along the line of sight, see 
Fig. [Hand Table H 

6.6 Evolution of flux densities 

The source of synchrotron luminosity is assumed here 
to be the flux of kinetic energy, L^, 



1 



2t/3 



P^ttR^V 



(94) 
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Radius (mas) 



Fig. 17 Cut of the mathematical intensity / of the rim model ( equation ((H31)) through the center (dotted line) of 



SN 1993 J and real data (empty stars). The parameters are a — 2 mas ,6 = 4 mas 



I(y=0)- 



1.73 and x = 125.19. 



T — ' — ' — ' — ' — r 





_J I I I I I t I > I I I > > I I I > I . I > I I > L. 

12 3*5 

Radius (mas) 



Fig. 19 Cross-section through the mathematical intensity / of ID diffusion in the optically thin case (dashed line, 
= 80.64) and real data (empty stars). Parameters as in Tabled 
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Table 3 Simulation of SN 1993J with 3D diffusion, opti- 
cally thick case with Ka =0.2. 



i{y=o) 
/(y=o) 



meaning 


value 


radius internal sphere 


2.01 (mas) 


radius of shock 


2.2(mas) 


radius external sphere 


5.0(mas) 


ratio observed intensities 


1.7926 


ratio optically thin case 


2.0491 


ratio optically thick case 


1.6741 




Fig. 18 Values of concentration computed with equa- 
tion H89p (full line) compared with the results of a Monte 
Carlo simulation (filled circles). The parameters in the 
Monte Carlo ID random walk are H^, = 0.5, Ai = 6, 
rn = 1.6510"^ pc and NTRIALS = 100. 



Table 4 Simulation of SN 1993J based on ID diffusion, 
optically thin case. 



symbol meaning 



value 



radius internal sphere 1.44(mas) 

radius of shock 2.79{mas) 

radius external sphere 4:.15{mas) 

ratio observed intensities 1.7926 

ratio optically thin 1.6430 



i{y=o} 



where R is the instantaneous radius of the SNR and 
p is the density in the advancing layer in which syn- 
chrotron emission takes place. The density in the ad- 
vancing layer is assumed to scale as R~'^, see formula 
([8]), which means that 



(95) 



The temporal and velocity evolutions can be given by 
the power law dependencies of equations p7)) and (PT|) 
and therefore 



Lrn (X t 



— ap(i+5ap — 3 



(96) 



The synchrotron luminosity L\ and the observed flux 
S\ at a given wavelength A are assumed to be propor- 
tional to the mechanical luminosity and therefore 



to 



(97) 



where 5*0 is the flux at t = to. The availability of 
synchrotron flux at 6 cm, see Table 1 and Fig. 14 in 
Marcaide et~d\ |2009l ). allows of fixing the parameters 
So- Fig. [20] reports the observed flux as a function of 
time as well as a theoretical evaluation using equation 
(|97l) . The astrophysical version of the above equation 



(98) 



5A=6crn =69.41 (-)-116 mjy , 

to 

with the time expressed in years. 



7 Conclusions 

Law of motion The first two parts of the law of mo- 
tion for SNR are thought to be a free expansion in 
which R (X t and an energy conserving phase in which 
R (X i^/s. A careful analysis of SN 1993 J in the first 10 
yr conversely suggest that R oc t^-^^. In other words the 
free expansion which follows the first Newton's law of 
motion does not corresponds to the observations. This 
observational evidence requires a law of motion which 
contains an adjustable parameter. Here we have cho- 
sen as a physical arguments the thin layer approxima- 
tion in a medium which has a decreasing density of the 
type oc R~'^ The two laws of motion here deduced are 
the classical nonlinear equation (jl3p and the complex 
relativistic- mechanics equation ((28)) . In the classical 
case of the thin layer approximation, we derived a use- 
ful asymptotic law of the type R oc t^/'^^~'^\ see equa- 
tion [TH This means that d = 1 — 1/(q:) and therefore 
the value of the d can be deduced from the observa- 
tional parameter a. Both the classical and relativistic 
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2 ' 4 ' 6 ' S 

time (years) 



Fig. 20 Observed time evolution of total flux densities of SN 1993J as given by VLBI data (empty stars) at 6 cm and 
theoretical relationship as given by formula (|97[) (full line), with d = 2.79 and ap = 0.82. The VLBI data at 6 cm are 
extracted from Table 1 in Marcaide et al. 2009. 



equations of the thin layer approximation can be used 
as a fitting function to deduce the remaining physical 
parameters , which are the initial radius and the initial 
velocity or /?. 

A third law of motion is deduced in the framework of 
the relativistic-hydrodynamics with pressure applying 
the momentum conservation, see equation p7p. The 
two relativistic cases here considered does not have an 
asymptotic solution. The evaluation of the merit func- 
tion I see equation as well the efficiency of the 
final radius, see equation ^ allows to eliminate the 
free expansion + the Sedov phase which are supposed to 
characterize the first two phases of the SNR's. The best 
results are obtained by the " ad hoc" piecewise function 
(|22p followed by the classical nonlinear equation ([T^]) 
once the ^ind efficiency are considered together. 

Relativistic velocities It is really necessary the 
relativistic treatment for the equation of motion? We 
briefly recall that at the moment of writing the max- 
imum observed velocity from Call H&K absorptions 
in SN2002bo is « 26000fcm/s as me asured at - 15 



days f rom B maximum , sec Figure 11 in lBenetti et < 
(j2004l ). A shift of 7 days characterizes the time delay 



between Gamma-ray burst ( GRB) and SN spectrum in 
the case of SN2003dh , see lMatheson et all (|2003[ ). A 
quadratic fit allows to extrapolate at t=-22 days from 
B maximum the velocity of expansion ; in particular 
we found sa 28710km/ s as a maximu m velocity for 
SN 20 02bo , see data in Figure lib of iBenetti et al. 
( 20041 ). This means that /3 ss 1/10.4 which is not far 



from the canonical P ~ 1/10 which marks the transition 
from classical to relativistic regime. 

Formation of the image The radial decrease in the 
density of the relativistic electrons from the position of 
the shock at 3 cm can be obtained by assuming a diffu- 
sive process and the resulting intensity of synchrotron 
emission can be calculated as the integral along the 



line of sight. Here, we considered the intensity profiles 
which arise from a 3D mathematical diffusion with con- 
stant diffusion coefficient in the optically thin case, see 
Fig. [131 In this process the particles can be acceler- 
ated by the Fermi I mechanism at the shock position 
in a region of thickness <C (c — 6) and {b — a). In this 
case the theoretical profile as given in Fig. [T3l toward 
the external region is concave up. A second model with 
diffusion in ID with a step length corresponding to the 
relativistic electron gyro-radius is also analyzed, see the 
concentration profile given in Fig. [15] and the intensity 
profile given in Fig. [TOl These diffusive processes allow 
to build up a theoretical 2D map of the shell-like inten- 
sity , see Fig. [TJ] A third model analyzes the radiation 
from a shell with constant density of emission, see Fig- 
ure 1171 This is the simplest model which produces an 
"U" profile in the cut of the intensity which toward the 
external region is concave down. In this case the ac- 
celerating mechanism can be the Fermi II mechanism 
characterized by multiple collisions in a shell having 
thickness b — a. 

A model based on the conversion into radiation of the 
flux of kinetic energy explains the observed decrease in 
flux at 6 cm, see equation[35]and Figure[2ni- This means 
that a direct conversion of the flux of kinetic energy 
which varies with time is a realistic model. Due to the 
fact that the observed proflle in intensity is concave 
up , see Figure [T3] , a direct acceleration through the 
Fermi I mechanism in thin region around the shock is 
an acceptable model. 
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